Re-Entrant Quantum Phase Transitions in Antiferromagnetic Spin-1 Ladders 
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In response to recent chemical attempts to construct higher-spin ladder materials from organic 
polyradicals, we study the ground-state properties of a wide class of antiferromagnetic spin-1 lad- 
ders. Employing various numerical tools, we reveal the rich phase diagram and correct a preceding 
nonlinear- sigma-model prediction. A variational analysis well interprets the phase competition with 
particular emphasis on the re-entrant phase boundary as a function of the rung interaction. 

PACS numbers: 75.10.Jm, 75.40. Mg, 75.40.Cx 



In 1983 Haldane Q awoke renewed interest in quan- 
tum spin chains predicting a striking contrast between 
integer- and half-odd-integer-spin Heisenberg antiferro- 
magnets. His argument was indeed verified in a spin-1 
material Ni(C2HgN2)2N02(C104) Q and was given an 
analytic support || as well. Since then the energy gaps 
in magnetic excitation spectra, that is, spin gaps, have 
been a central issue in materials science. In the last 
decade more and more researchers made a wide variety of 
explorations into the spin-gap problem, such as the spin- 
Peierls transition in inorganic compounds Q, quantized 
plateaux in magnetization curves q|, and antiferromag- 
netic gaps in the ferromagnetic background Among 
others Dagotto et al. fit] pointed out that another mecha- 
nism of the gap formation should lie in a ladder— two cou- 
pled chains. A spin gap was indeed observed in a typical 
two-leg ladder material SrCu203 ||. Moreover, super- 
conductivity was brought about in its hole-doped version 
(SrCa)i4Cu2404i [Bj. Ladder systems caused us further 
surprise exhibiting excitation spectra varying with the 
number of their legs [0 . 

So far metal oxides have been representative of ladder 
materials. Though molecule-based ones |ll],|l2| have been 
synthesized in an attempt to reduce the spin gaps and 
obtain experimental access to them, the situation of cop- 
per ions supplying the relevant spins remains unchanged. 
Therefore they are all spin-i antiferromagnets. In such 
circumstances there has occurred a brand-new idea of 
constructing purely organic ladder systems. Katoh et al. 
p3[ synthesized novel organic biradicals and tetraradicals 
which crystallize to form an antiferromagnetic ladder of 
spin-i and that of effectively spin-1, respectively. Their 
polyradical strategy has yielded further harvest such as 
an effective spin-1 antiferromagnet on a honeycomb lat- 
ticeJ14j and a ladder ferrimagnet of mixed spins 1 and 
i |15|], displaying the wide tunability of the crystalline 
structures in higher-spin systems as well. 

Distinct spin-gap mechanisms may lie in higher-spin 



ladders and quantum competition between them must 
lead us to further enthusiasm for ladder systems. There 
exist pioneering works in this context. Sierra Jl6| gen- 
eralized the well-known nonlinear-sigma-model analysis 
fli"|,[i~7t on quantum spin chains to multi-leg ladder sys- 
tems. His findings suggest that only an odd number of 
half-odd-integer-spin chains are massless, supporting ex- 
perimental observations on a series of spin-^ ladder an- 
tiferromagnets [p|JT8| . The technique was further devel- 
oped for spatially inhomogeneous ladders JTq| . Mixed- 
spin ladders [ p0pl| were also investigated with partic- 
ular emphasis on the competition between massive and 
massless phases. 

In comparison with extensive calculations j7 22-24|] 
on spin-i ladders, there exist few quantitative investi- 
gations into higher-spin ladders [g5| and most of the 
above-mentioned predictive theories remain to be veri- 
fied. Hence, in this article, we solve the ground-state 
properties of antiferromagnetic spin-1 ladders with two 
legs. Employing various numerical tools and comple- 
menting them by a variational argument, we elucidate 
the valence-bond-solid-like nature of their ground states, 
which is in contrast with the spin-liquid or rcsonating- 
valence-bond ground states |2^]23|,^6|,^7j of spin-i two- 
leg ladders. The obtained phase diagram is reminiscent 
of the preceding sigma-model prediction |l^] but con- 
tains a re-entrant phase boundary, which can never be 
extracted from any field-theoretical argument. 

Considering that an advantage of assembling organic 
open-shell molecules into a magnetic material is the 
isotropic intcrmolecular exchange couplings, while the 
polyradical strategy is accompanied by spatial variations 
in magnetic interaction fl4|| , we treat a wide class of spin- 
1 antiferromagnetic ladders 
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FIG. 1. Phase diagrams for the antiferromagnetic spin-1 
ladder with two out-of-phase legs, (a) A field-theoretical pre- 
diction [19] . The two critical lines (dashed lines) derived from 
the effective sigma model for ladders are inconsistent with 
the sigma-model analysis on isolated chains (•). They re- 
main far apart from each other even in the decoupled-chain 
limit r = 0. Therefore qualitatively patched-up phase bound- 
aries (solid lines) were predicted, (b) Our numerical findings. 
The series-expansion estimates are shown by x, while the 
level-spectroscopy analyses by □ (L = 6) and o (L = 8). 
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FIG. 2. Demonstration of the level spectroscopy. The low- 
est-lying two eigenvalues in the subspace of zero magnetiza- 
tion as functions of r cross at transition points provided the 
twisted boundary condition is imposed on the Hamiltonian. 

where the bond-alternation parameter 7jj is defined in 
two ways as 

_ f 1 + {-l) %+j 5 OP (out-of-phase legs) , , , 
~ \ 1 + {-iy8 w (in-phase legs) . [ > 

We calculate the region of < 5op(fa>) < 1 and here- 
after set J\/J\\ to r(> 0). Martin-Delgado, Shankar, 
and Sierra |1!| studied the cases of out-of-phase legs de- 
riving a low-energy-relevant sigma model. For the spin-S* 
ladders with two out-of-phase legs, the topological angle 
in the effective sigma model turns out 8irS5op/{r + 2) 
and reads as the critical lines 8S5qp — (2n + l)(r + 2) 



(n — 0, ±1, • ■ •). However, these findings do not smoothly 
merge with the well-established critical behavior in one 
dimension, 25(1 — 5) = 2n + 1 |Q, as is shown in Fig. 
|l|(a). Thus, it is necessary to verify the true scenario all 
the more in higher dimensions. 

One of the most reliable solution may be a numerical 
analysis p8| on the phenomenological renormalization- 
group equation f29|j . However, the scaled gaps are ill- 
natured due to the close critical points, so as to make 
the fixed points hard to extract from available numer- 
ical data. Then we switch our strategy to the level 
spectroscopy [Q, the core idea of which is summarized 
as detecting transition points by crossing of two rel- 
evant energy levels. Although the method is generi- 
cally applicable to the Gaussian critical points [plf , no 
explicit change of symmetry accompanies the present 
phase transitions and therefore any levels do not cross 
naively. In order to overcome the difficulty of this 
kind, Kitazawa [|32"f proposed the idea of applying the 
twisted boundary condition, that is to say in the present 
case, setting the boundary exchange couplings equal to 

- Eli ^llT^^^i+^^i-^^i)- Thenthe en " 
ergy structure of the Hamiltonian is changed and the low- 
est two levels are led to cross at transition points, which 
is demonstrated in Fig. ||. Due to the limit of time and 
memory well spent, we have restricted our calculations up 
to L = 8. We plot in Fig. 0(b) bare findings for the cross- 
ing points at L = 6 and L = 8 rather than extrapolate 
them trickily. We are sure that the data uncertainty still 
left is within the symbol size. A series-expansion tech- 
nique |3^j3^] guarantees the level spectroscopy to work 
well. Starting with decoupled singlet dimers on legs or 
rungs and expanding the energy gap as a power series 
in a relevant perturbation parameter, we can obtain a 
partial knowledge of phase transitions. Here we have 
calculated the gap up to the ninth order and further ap- 
plied the Dlog Pade approximants [^5| to them. The 
thus-obtained phase boundaries, which are also shown 
in Fig. 0(b), elucidate the nature of the phase compe- 
tition, that is, the Afffeck-Kennedy-Lieb-Tasaki (AKLT) 
valence-bond-solid (VBS) || on a snakelike path versus 
decoupled dimers. 

The most impressive findings are re-entrant quantum 
phase transitions with increasing r. The preceding sigma- 
model analysis Jl9| is indeed enlightening but never able 
to reveal this novel quantum behavior. In order to char- 
acterize each phase, let us consider a variational ap- 
proach. We know that singlet dimers on rungs [Fig. ||(h)] 
are stabilized for r — > oo, whereas either dimers on legs 
[Figs. |(d) and |(e)] or the double AKLT VBS [Fig. |(c)] 
for r — > 0. Two more interchain VBS states [Figs. J5](i) 
and ^(g)] may be adopted as variational components for 
the intermediate-r region. Thus the linear combination 
of Figs. ||(c) to ||(h) can be an approximate ground-state 
wave function for spin-1 ladders. Since the present vari- 
ational components are all asymptotically orthogonal to 
each other, the variational ground state turns out any 
of them itself |3(|. The thus-obtained phase diagram is 
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presented in Fig. ||. The significant stabilization of the 
intermediate phase, which is now characterized as SH, 
and the resultant re-entrant phase boundary are success- 
fully reproduced. Considering that a couple of critical 
chains immediately turn massive with their rung inter- 
action switched on [Q, the point C should coincide with 
the point A under more refined (and thus inevitably nu- 
merical) variational investigation. 

The present variational calculation implies possible 
phase transitions for in-phase-leg ladders as well, but this 
is totally due to the naive wave function. Numerical ob- 
servation of the energy structure ends up with no gapless 
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FIG. 3. Plaquette-singlet-solid and valence-bond-solid 
states relevant to the two- leg antiferromagnetic spin-1 lad- 
ders. • denotes a spin | and their segment linkage means a 
singlet formation. O represents an operation of constructing 
a spin 1 by symmetrizing the two spin ~'s inside. 
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FIG. 4. Variational phase diagrams for the two-leg an- 
tiferromagnetic spin-1 ladders. The thick solid lines de- 
scribe phase transitions, whereas the thin ones represent the 
crossover of the ground-state nature within the present vari- 
ational scheme. The dotted line is only a guide for eyes. 



point in this region. The sigma-model approach also 
concludes no critical point, giving the topological an- 
gle AnS independent of both r and 5p . The key to the 
ground-state nature of in-phase-leg ladders is the four- 
spin correlation pf|. Let us consider interacting four 
spins of S — ^ which are described by the Hamiltonian 

H = Yh=\ Si,r Si,i+i+rJ2 l jt?i Si.j-S 2 .j ■ In terms of the 
Schwinger boson representation: S + = a^b; S~ = atf; 
S z = |(a T a — tfb); S = ^{a^a + b^b), their ground state 
is explicitly given as 



Mr)) i = 



cos0(r)n(«L^ + i- fo l«M + i^ 

|0>, , (3) 



i+i 



3=1 



where |0); is the Bose vacuum and 9(r) is given by 
tan 9(r) = r — 1 + \/l — r + r 2 . As r varies from to oo, 
9 moves from to ? , that is, the leg dimers continuously 
turn into the rung dimers. Using the plaquette singlet 
state (^), we can in principle construct much better vari- 
ational wave functions for spin-1 ladders particularly in 
the in-phase-leg region as 



2/-1 



x cos (p(Sip)\ip{r)) 21-1 +sin^(5ip)|^(r))2i 



(4) 



where Vij represents an operation of symmetrizing the 
two spin i's at site (i, j) into a spin 1. Now that 9(r) 
is a continuous function of r and may here deviate from 
that in Eq. (|^) , a naive optimization |36| of Eq. (|4|) is no 
more feasible. However, the refined variational scheme 
shows us more. For better understanding of the wave 
function (|J), we visualize in Fig. |^ its special forms for 
4> = 5 (a) and <fi — (b), which are most stabilized at 
<5ip = and Sip — 1, respectively. The snapshots of PSS 
(DSP) at 6 = 0, 9 = f , and 9 = § are nothing but 
the variational components DH (IPLD), DP, and RD, 
respectively. Through the plaquette singlet resonance, 
any snapshot of PSS (DPS) can turn into another with- 
out any explicit transition. PSS and DPS share DP and 
RD as their snapshots. Thus, the extended variational 
wave function (Q) is expected to erase the artificial first- 
order transition lines (thin solid lines), reducing the only 
discontinuity wall AB to a point. The point A belongs 
to the same universality class as the spin-i Heisenberg 
chain j3l] . We are all the more convinced of the imme- 
diate gap formation with r moving away from 0. On the 
other hand, neither PSS nor DPS includes both OPLD 
and SH and therefore the two critical lines in the out- 
of-phase-leg region survive against the plaquette singlet 
formation. 

The generalized string order parameter |37j 0{9) = 
lirri| i _ J |^ 00 (Sf YljZi exp[i# 5f ]Sj) is also useful in char- 
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acterizing the ground state. 0(9) distinguishes between 
one-dimensional VBS states by its 9 dependence @. 
Hence, measuring it on the linear-chain and snake paths, 
we can detect the transitions between DH, OPLD, SH, 
and RD, as is shown in Fig. 0. If we specify the 
transition through a change of the 9 dependence in the 
vicinity of 9 = 7r, that is, the change from the con- 
vex curve to the concave one, we obtain the transi- 
tion points (S OP ,r) = (0.245,0) and (0.6,2.125), which 
are in good agreement with the numerical findings in 
Fig. |l|(b). Whatever path we take for 0(9), its peak 
never sits on 9 = tt in the in-phase-leg region, sug- 
gesting that we can not observe the Haldane state of 
any kind there. The plaquette singlet formation can 
instead be visualized by extending 0(9) to ladders as 
lim| i _ iHoo (^ )i 5| ji FID 1 eM^S^+Sl^SIJ |§. 
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FIG. 5. Quantum Monte Carlo estimates of the generalized 
string order parameter defined on the linear-chain (a) and 
snake (b) paths. The dashed lines represent the analytic cal- 
culations for the AKLT VBS (| sin 2 |) and decoupled dimers 
(| sin 2 9). Sop runs on the line of r = 0.0 in (a), while r runs 
on the line of <5op = 0.6 in (b). 

We summarize our rich harvest as the plaquette sin- 
glet formation and the re-entrant quantum phase tran- 
sition via the snake Haldane state in staggered spin-1 
ladders. The re-entrant phase boundary is peculiar to 
spin-1 or possibly integer-spin ladders. Half-odd-integer- 
spin chains are generically critical at the translationally 
symmetric point 6 = 0. Therefore, the critical line can 
not exhibit the initial re-entrant behavior in the r-S phase 
diagram, because it should be symmetric for ±6. Chem- 
ical tuning of bond-alternating critical chain compounds 
p8[ , as well as the organic-radical-based materials re- 
search is strongly encouraged. 
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